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ABSTRACT 



We apply the theory of parton-parton total cross sections at large s, due to 
Lipatov and collaborators, to compute the inclusive cross section for jets which 
accompany a large s parton scattering process. 
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The enormous dynamic range for QCD processes opened by the Tevatron col- 
lider has given a new impetus to detailed investigations of the dynamics of quark 
and gluon scattering. Though the most prominent experimental investigations 
have involved the cross section for 2-jet production, there are also new results on 
3- and multi-jet processes. In these higher-order processes, there are two or more 
different Lorentz invariants which set the scale of the momentum transfer, and, 
typically, these invariants may differ by large factors. The study of QCD processes 
with large ratios of invariants — semihard processes — is complicated theoretically 
because it typically involves the resummation of infinite classes of Feynman di- 
agrams. It is complicated experimentally because it requires jet detection in a 
large rapidity interval. Nevertheless, this regime is a fascinating one, with new 
and nontrivial applications of QCD perturbation theory. In addition, as Bjorken 
has recently emphasized, the study of jet correlations over large rapidity intervals 
may offer interesting signatures of new physics! 11 

In this paper, we will discuss the following situation, which involves jet dynam- 
ics with a large ratio of invariants. Consider a process in which gluons or quarks 
scatter with large center-of-mass energy and only moderately large transverse mo- 
mentum. Call the squared center-of-mass energy of the original parton-parton 
scattering s, and let m be a typical value of the final parton transverse momenta. 
In the semihard regime, s » m 2 » Aq CD , Lipatov and collaborators' 2 41 have 
shown that, in this regime, the rapidity interval between the scattered partons is 
filled in by the radiation of additional gluons, roughly uniformly spaced in rapidity, 
all with transverse momenta of order m. In principle, the properties of these radi- 
ated gluons can be observed by choosing events with two jets at large positive and 
negative rapidity and then measuring the production rate for jets at intermediate 
rapidities which accompany these jets. In this paper, we will call the jets at large 
rapidity tagging jets, and we will call the jets that accompany them minijets. A 
few results on the spectrum of these accompanying gluon jets are known from the 
work of Balitsky, Fadin, Kuraev, and Lipatov (BFKL), refs. 2-4, and from suc- 
ceeding work of Levin and Ryskin. M In this paper, we return to this question and 
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give a systematic procedure for computing the inclusive spectrum of accompanying 
jets, along with numerical estimates for Tevatron and SSC energies. 

Total Cross Section for Tagging Jets 

In order to understand the spectrum of minijets which accompany a set of 
tagging jets, we should first review the QCD prediction for the total cross section 
for producing these tagging jets. This prediction has been presented by Mueller and 
Naveletf 1 using the expression for the asymptotic parton-parton total cross section 
which is the principal result of the BFKL analysis. We will use the ingredients in 
this prediction to construct the related inclusive cross sections. 

Following Ref. 6, we consider the scattering of two hadrons of momenta pa 
and pb in the center-of-mass frame, with the z axis along the beam momenta, and 
we imagine that we tag two jets at the extremes of the Lego plot, with the rapidity 
interval between them filled with minijets. The tagging jets can be characterized by 
their transverse momenta and by their longitudinal fractions xa, %b with respect 
to their parent hadrons. It is simplest to consider the cross section for producing 
two tagging jets with transverse momenta greater than a minimum value m. Then 
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(1) 

crtot(s), 



where s = 2pA -pb^A^b is the parton-parton squared center-of-mass energy, and 
otot is the BFKL total cross section for gluon-gluon scattering, which we will dis- 
cuss in a moment. Eq. (1) includes the effects of quarks using the observation of 
Combridge and Maxwell that, in a process with large rapidity intervals, the leading 
contribution to any scattering process comes from gluon exchange in the crossed 
channel. In writing (1), we assume that the values of xa, %b are sufficiently large 
that the parton distributions can be computed from ordinary Altarelli-Parisi evo- 
lution; semi-hard QCD adds additional complications when these fractions become 
small! 81 
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The core of eq. (1) is the BFKL function crtot(s), which is given by 



crtot(s) 



8 nN 2 a 2 s 
N 2 -l 2m 2 



F(Y), 



(2) 



where N c = 3 is the number of colors in QCD, Y = log(s/m 2 ), and F(Y) is a 
dimensionless function which we will discuss below. The strong coupling constant 
is evaluated at a scale m 2 ; the running of a s is sub leading in the BFKL theory. 
Setting F — 1 gives the cross section for gluon-gluon scattering at the lowest order 



The BFKL theory systematically corrects the lowest-order QCD result for a tat 
by summing the leading logarithms of (s/m 2 ). This is done in three stages, as 
shown in Fig. 1. First, one simplifies the lowest-order QCD diagrams for multigluon 
production, shown in Fig. 1(a), for the case in which the emitted gluons are widely 
separated in rapidity. The gluon emission vertex is replaced by a non-local gauge- 
invariant effective vertex' 21 . Next, one sums the leading corrections to the forward 
amplitude with color octet in the t-channel, as shown in Fig. 1(b). The result has 

roi 

the form of a Regge pole with an infrared- sensitive trajectory. Finally, one uses 
this resummed, effective gluon exchange to compute the forward amplitude with 
color singlet in the t-channel! 41 The imaginary part of this object is the gluon-gluon 
total cross section. 

The last two step of this procedure involve the solution of integral equations 
constructed by BFKL. In the final step, the integral equation can be solved ex- 
plicitly for the imaginary part of the forward amplitude to find tagging jets with 
transverse momenta k^±, ks±, both greater than m, separated by a rapidity in- 
terval Y. The solution of the equation involves a Laplace transform with respect 
to rapidity. Then the amplitude is given by 



of QCD, integrated over transverse momenta with \ki±\ > m, in the limit s >> m 2 . 



a+ioo 




(3) 
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where the Laplace transform has the representation 

i r (u2 \-l/2+iv(u2 \-l/2-iv 

(2t) „ J u-u(n,u) 

In this equation, 0,4 — 0b is the azimuthal angle between the transverse momenta 
of the tagged jets, and 

u(n, v) = ty(l) - Re^(^±i + ii/)] , (5) 

7T Z 

with -0(z) the standard logarithmic derivative of the Gamma function. Throughout 
this paper, we ignore azimuthal correlations and keep only the leading, n = 0, term 
of (4). Near v — 0, uo(v) = uj(0, v) has the expansion 

u (v) = A-Bu 2 + ---, (6) 

with 

A = ^*Uog2, B=l^iC(3) (7) 

7T 7T 

The integral over a; in (3) can be done easily by picking up the pole. Then we 
may integrate over the jet transverse momenta to obtain the enhancement factor F 
in eq. (2). The integrals over k^± and ks± are singular and depend on the cutoff 
m; this gives the factor to~ 2 in (2). Comparing with that formula more closely, we 
find 

00 

— oo 

The exponential growth of F(Y) with the rapidity interval is associated with mini- 
jet production. Using (6) to expand about the saddle point at v — 0, we can see 
that F(Y) has the asymptotic behavior 

F(Y) ~ T with z = -^Y. (9) 

V7C(3)7r*/2' 

Mueller and Navelet showed that this asymptotic form is an accurate representation 
for z > 0.2. 
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Mini-jet inclusive cross section 

The BFKL total cross section is a sum over multi-jet emission processes. Thus, 
it is not difficult to pull this amplitude apart and find the contributions from final 
states with a jet in a specific region of phase space. The cross section for producing 
a jet at rapidity y and transverse momentum gj_ in association with the tagging 
jets described in the previous section may be computed as the sum of diagrams 
shown in Fig. 2. We use the BFKL propagator (3) to represent the ladder and 
Lipatov's vertex from ref. 2 to represent the gluon emission. For the moment, we 
retain the dependence on the transverse momenta kA±, k B ± of the tagging jets. 
Then the inclusive cross section for minijet emission is given by 



1 da 



a to t dyd 2 q_ L d 2 k A ±_d 2 k B1 _ 



16a s N c m 2 f d 2 k 1± d 2 k 2 ± . , 2x ( 2 ) n , \ 
: n FiY) J W W? ( } ( 1X " 2± ~ ^ 

■ , 2 \ 2 f( k A±,ki±,yA-y)f(k 2± ,k B± ,y-y B ), 

(10) 

where y A , y B are the rapidities of the tagging jets, with Y = y A — y B . This formula 
has been derived previously by Levin and Ryskin. lJ ' 

As before, we ignore all correlations in azimuthal angle by integrating or av- 
eraging each transverse momentum over its angle 0. Note that the correlations 
in are, in any event, subleading. With this simplification, we can represent the 
momentum delta function in (10) as 



(27r) 2 5 ( - 2 \k l± -k 2± -q ± ) = J d 2 be Mk^-k 2 ^) = J d 2 b J ( bkl± )J (bk 2± )Jo(bq ± ). 

(11) 

By introducing the representation (4), we can now perform the integrals over the 
four k{±. Notice that the integrals over k A j_, k B ± depend on the infrared cutoff, 
as in the derivation of (8). However, the integrals over ki±, k 2 ± converge in the 
infrared. The resulting expression is a function of two Laplace transform variables, 
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which we will call vi, v 2 . We are left with the integral over 6, which reduces to 

oo 

J dbb- 1+2i ^-^ J (b qi _) (12) 



To make this integral well-defined, we deform the contours of integration over the 
Vi so that lm(u 2 ) > lm(u\). We must keep |Im(i/j)| < \ in order that the integrals 
over the ki± converge. We will show below that there is a saddle point for the V{ 
with this property. 

After performing the integral over b, we are left with the following expression 
for the inclusive cross section for the production of one minijet with rapidity y and 
transverse momentum q in association with tagging jets with transverse momentum 
greater than m: 



oo 

I da _ 1 r dv x dv 2 u { v l){yA - y) M(y- yB ) N cU S( m 2 , t(«/i-m)/ -i \ 
a tot dyd 2 q ± F(Y) J 2tt 2tt n 2 q 2 V } \u 1 -u 2 ) 



-oo 



l 2 _ 

1 VQ 1 -iu 1 )Y{l + i(v l -v 2 ))Y{\ + iu2) 1 

;i - ivr) r(i + iv x ) r(i - i{v x - u 2 )) v{\ - iv 2 ) {\ + %v 2 ) 

(13) 

The variables V{ are integrated along contours with ^ > lm(u 2 ) > lm(ui) > —\. 

The same method yields an integral expression for the n-minijet inclusive cross 
section. The integration variables are the n + 1 Laplace transform variables of 
the intermediate BFKL propagators. By the argument given above, these are 
ordered along the imaginary axis from the tagging jet B to the tagging jet A: 
\ > Im(z/ n+ i) > ■ ■ • > lm(ui) > —\. The full inclusive cross section is given by 

1 da 



°tot n?=i d2( li± d yi 

1 n + l r a,.. n ( AT „, ,„2 



F(Y) 11 J 2vr 11 1 n\ 2 ± Km? ) 

1 Y{\ - ivi) / A T[l + % {vj - n+i)] \ V{\ + ivn+i) 1 

11 rn - 



(i - r(i + %v x ) \ f \ r[i - %{vi - v i+1 )\ / r(i - w n+1 ) (± + w n+1 ). 

(14) 
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Asymptotic Evaluation 

In eq. (9), we simplified the integral formula for the total jet-jet cross section 
by evaluating the integral asymptotically around an appropriate saddle point. It 
is straightforward to evaluate the one-jet inclusive cross section asymptotically in 
the same way. We will find a joint saddle point in the variables vi, v 2 and expand 
the double integral around this point. The expansion will be valid in the limit in 
which (yA — y) and {y — Vb) are both large. We neglect the terms in the square 
brackets in eq. (14), since at asymptotic energies the saddle point will be very 
close to the origin in the plane v\, v 2 . 

The joint saddle point will occur with both of the V{ on the imaginary axis: 
V{ = —irii. Using the expansion (6), we find the saddle point conditions: 



2(jM - y)Bni - log(gi /m 2 ) - 
2(y - y B )Bn 2 + \og{q 2 Jm 2 ) + 



1 



n\ - n 2 
1 

ni - n 2 







(15) 



From these equations, 



ni = (y- ubW, n 2 = -(va - y)v, 



(16) 



with 



and 



V 



71 \og(q 2 L /m 2 ) 



5QN c ((3)a s (y A - y)(y - vb) 



[1 + a], 



a = Jl + 



ll2N c ((3)a s (y A - y)(y - y B ) 



7TY\n 2 (ql/m 2 ) 
Expanding about this saddle point, we find 



(17) 



(18) 



da 



(Jtot dyd In q]_ 
with a as in (18). 



N c a s 1 
^T~a eXP 



a-1 



+ 1/2 



(19) 
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According to this equation, the jet inclusive cross section, falls off faster than 
the scale invariant dependence (d 2 q±/q^_) expected for an approximately fixed value 
of a s . At small transverse momentum, where we can ignore the log(g^_/m 2 ) terms 
in solving (15), there only is a small modification by a factor (q^)~ Yri . At large 
transverse momentum, the cross section falls off faster than any power of q±. This 
latter, doubly asymptotic, limit has been found previously by Ryskin using another 
method. 

Numerical Evaluation 

The formula (13) is also an appropriate starting point for a numerical eval- 
uation of the one jet inclusive cross section. It is important to study the exact 
behavior of this formula numerically, for two reasons. First, the total rapidity in- 
terval for jet-jet scattering processes is not large at energies now available, and so 
it might be troublesome to divide this interval into several pieces, each of which 
must be large to justify an asymptotic analysis. Second, this problem is exacer- 
bated by the presence of the poles in (13) which restrict the region of the imaginary 
axis through which the v\ contours can pass. At asymptotic energies, the saddle 
points found in the previous section lie very close to v\ — 0, and so these poles are 
not relevant. However, the evaluation of the saddle point locations with realistic 
parameters puts n\ and ri2 close to Thus, we might expect large discrepancies 
for the saddle point result in realistic cases, and, indeed, we will find this below. 

To compute the minijet cross section explicitly, we evaluated the double integral 
in (13) numerically using the contours Imz/i,Imz/2 = ±|. We evaluated a s = 
a s(<lA-), scaled from a s {mz) = 0.12 using 1-loop evolution with 5 flavors. 

For (8), we saw that the saddle point approximation was an accurate one. 
However, this turns out not to be true for the asymptotic evaluation of (13). In 
Fig. 3, we compare the exact and asymptotic formulae for the one jet inclusive 
cross section for gluon-gluon scattering processes at center of mass energies from 
1 to 10 9 TeV, with m = 20 GeV. There are large discrepancies between the two 
formulae which disappear only extremely slowly, as the inverse of the logarithm of 
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the energy. In fact, the asymptotic formula is never accurate. Since the total cross 
section (9) rises as a power of s, the contribution from the exchange of one gluon 
ladder must at some point be unitarized by contributions from multiple ladder 
exchanges. This correction is already important at 100 TeV, and changes the one 
jet inclusive cross section essentially at higher energies. 

However, even if the asymptotic formula is not valid, we can estimate the 
minijet inclusive cross section directly from (13). It is important to recall that 
this equation, as well as (8), involve additional asymptotic approximations. The 
BFKL integral equation is derived under the assumption that the emitted gluons 
are widely separated in rapidity. The solution of the color octet exchange problem 
shown in Fig. 1(b) involves the assumption that the leading Regge pole in the u 
plane dominates over the remaining singularities. However, there is no reason why 
either of these effects should lead to large corrections under realistic conditions. 

With this set of approximations, then, we present in Fig. 4 predictions from 
(13) for the minijet inclusive cross section at Tevatron and SSC energies, taking 
the longitudinal fractions xa, ib of the tagging jets momenta to be 0.1. The 
cross section is computed for an observed jet in the center of the rapidity interval 
between the tagging jets. However, the minijet cross section is almost independent 
of the position of the jet in the rapidity interval. These curves, and also those of 
Fig. 3, show a shift of the transverse momentum distribution toward higher values 
as the energy of the original scattering process — and, therefore, the length of the 
gluon cascade — increases. This effect can also be seen in eq. (19), although the 
dependence of that formula on q± and y does not match our quantitative results. 

The final situation in which we are left is somewhat ambiguous. On one hand, 
we have improved the understanding of the physics of the BFKL total cross section. 
We have presented predictions from the BFKL theory for the one minijet inclusive 
cross section, and we have presented a formula which can be used to evaluate 
higher order jet correlations. However, we have found that the natural asymptotic 
evaluations of these formulae are not accurate at the energies of present colliders. 
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It will be interesting to see whether the estimates presented in Fig. 4 will be 
confirmed experimentally. However, we did not succeed in presenting a quantitative 
and characteristic prediction which can be used to confirm the BFKL resummation. 
That remains an interesting problem for the future. 
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FIGURE CAPTIONS 



1) BFKL resummation: (a) emission of gluons; (b) construction of the octet 
amplitude; (c) construction of the singlet amplitude. 

2) Sum of graphs leading to the inclusive jet cross section. 

3) Comparison of the exact and asymptotic evaluations of the one jet inclusive 
cross section as a function of the jet transverse momentum, at the parton 
center-of-mass energies y/s = 1, 10 3 , 10 6 , 10 9 TeV, and m = 20 GeV. The 
upper curves represent the exact evaluation of integral (13), and the lower 
curves its asymptotic evaluation (19). 

4) One jet inclusive cross section as a function of the transverse momentum 
of the jet for typical Tevatron and SSC energies. We choose y/s = 0.18 
and 4 TeV respectively for the Tevatron and the SSC parton center-of-mass 
energies, and we label the curves according to the minimum jet transverse 
momentum. 
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